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ABSTRACT 

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital 
decay and mass evolution of massive black hole (MBH) pairs down to scales of ~ 30 pc during minor mergers 
of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback 
from both processes. We consider 1:10 merger events starting at z ~ 3, with MBH masses in the sensitivity 
window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics 
and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive 
MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes 
of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside 
the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time 
the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of 
the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become 
as large as 1 :2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair 
itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down 
to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not 
necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay 
between gas accretion and merger dynamics. 

Subject headings: black hole physics — cosmology: theory — galaxies: evolution — hydrodynamics — meth- 
ods: numerical 



1. INTRODUCTION 



The ubiquity of massive black holes (MBHs) a t 

i et al.1 [T998b together 



the cen- 
ters of galactic spheroids dRichstone t 
with the "bottom-up" nature of galaxy formation in the cur- 
rently favored ACDM cosmology (e.g.. lWhite & Ree si [19781) 
suggest that MBH pairs at sub-kpc scales should form in 
galactic nuclei during the hierarchical assembl y of structures 
dBegelman et al.lfl980h see the recent review by Colpi & Dottil 
2009). Such pairs may undergo further orbital decay, become 
Keplerian binaries, and eventually coale sce via the emission 
of gravitational waves (Haehnelt 1994), which will be one 
of the main targets of the next generation of gravitational 
wave detector s such as the L aser Interferometer Space An- 
tenna (LISA) dVecchiol 120 04). Moreover, observations indi- 
cate that the masses of MBHs correlate with properties of their 
host spheroids, i ncluding the luminosity , mass, and velocity 
dispe r sion (e.g.. |Ma gorrian & al. 1998; Ferrare se & Merrita 
120001: iGebhardt & al N2000I) . Such relations suggest the ex- 
istence of fundamental physical mechanisms that link MBH 
assembly and galaxy formation, and may connect the proper- 
ties of galaxy mergers with the resulting MBH binaries. 
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Given all these facts, the study of the dynamics of 
MBHs during galaxy mergers becomes especially impor- 
tant as a means to connect the cosmological assembly of 
galaxies with that of MBH pairs and binaries. Many nu- 
merical studies have focused on the connection between 
galaxy mergers and MBH growth, specifically in relation 
to the MBH final ma ss, scaling relations and phenomenol- 
ogy o f quasars (e.g.. iDi Matteo et al. 2005; Hopkin s et alj 



2005; Boylan-Kolchin_&Ma| 120071; lYounger et al.1 120081; 
Johansson et al. I2009h . However, much less attention has 
been devoted to investigating the orbital decay and evolution 
of both MBHs during mergers, particularly in the unequal- 
mass re gime which comprises the vast majority of such 
events, kazantzidis eTaT] d2005h and ICallegari et all (|2009) 
(hereafter C09) showed that the formation of unequal-mass 
MBH close pairs at sub-kpc scales is sensitive to the details 
of the gas dynamics during the merger process. In particular, 
they found that, due to the combination of gas dynamics and 
star formation, a pair of MBHs can form efficiently in 1:10 
minor mergers at scales below 100 pc, provided that galaxies 
are relatively gas-rich and that the mergers occur at relatively 
high redshift (z ~ 3) when halo densities are higher and dy- 
namical f riction timescales co rrespondingly shorter, accord- 
ing to the IChandrasek har ( 1943) standard formulation. How- 
ever, these papers did not model the accretion onto the MBHs, 
the evolution of their mass ratio and its dependence on the 
merger dynamics. Such investigations are important, as the 
mass ratio of MBH pairs is a fundamental parameter which 
drives their e volution on nuclear scales (e.g., Dot ti et alj2007l ; 
lLodato et al.ll2009l) . 

In this paper, we explore in detail for the first time the for- 
mation and mass evolution of close MBH pairs down to scales 
of ~ 30 pc using controlled smoothed particle hydrodynam- 
ics (SPH) simulations of minor mergers between disk galax- 
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ies. The merging systems have a mass ratio of 1 : 10 and 
our simulation suite includes the effects of star formation and 
accretion onto the MBHs, as well as feedback from both pro- 
cesses. We model merger events associated with MBH pairs 
whose gravitational wave emission should be in the sensitiv- 
ity window of LISA: they have masses ~ 10 5 M Q around the 
epoch of a predicted peak of MBH pair formation in merger s 
happening at z ~ 3 (IVolonteri et al.ll2003l:ISe"sana et al.ll2005l) . 

2. NUMERICAL SIMULATIONS 

Our reference galaxy model is a z = 3 progenitor of a Milky 
Way-type disk galaxy, with a virial velocity y v ; r = 145 km s" 1 . 
Such a progenitor is c onstructed as suming a constant VW be- 
tween redshift and 3 dLi et al.ll2007l) . and rescaling the mass 
and size of a z = model by the ratio [H(z = 3)/H ] = 0.224 
of the Hubble constant at z = 3 over its present-day value 
dMo etal.ll 19981) for a ACDM concordance cosmology (Hq = 
70 km s" 1 Mpc -1 , O m = 0.3, fl\ = 0.7). The model consists 
of three components: a dark matter halo, a stell ar and gaseous 
disk, a nd a stellar bulge. The halo is a spherical Navar ro et alJ 
( 1996) model with a virial mass of Myj r = 2 3 x 1O H M , a spin 
parameter A = 0.04 {Vitvitska et al. 2003), an d a concentra- 
tion c = 3 consistent with the relations found by lBullock et alJ 
d2001l) . The halo has been adiabatic ally contracted to respon d 
to the growth of the disk and bulge (Blumenthal et alj|1986h . 
The baryonic disk follows an exponential profile, with a to- 
tal mass Md = 0.04M v i r , a radial scale-le ngth R c { = 1.1 kpc 
determined according to Mo et alJ (1 19981) . a scale-height of 
Zd = O.IRd, and a m ass fraction in gas denoted by /„. The 
bulge is a spherical Hernquist (1990) model with a mass 
Mb = 0.008M v ; r , yielding a ratio of bulge-to-total stellar mass 
B/T ~ 0.2 (depending on / g ), and a scale radius a& = 0.2/^. 
The satellite galaxies are constructed following the same rela- 
tions between structural parameters and a mass in each com- 
ponent scaled down by g ga i = 0.1, with a disk scale-length of 
300 pc. 

The A^-body realizations consist of 10 6 particles in the dark 
matter halo and 2 x 10 5 stellar particles in each galaxy. In 
addition, each galaxy is initialized with 10 5 gas particles in 
the disk, except for one case (see Table Q}. We adopted a 
force resolution of 45 pc for both the dark matter and baryonic 
particles of the larger galaxy, while for the satellite galaxy we 
used 20 pc. 

A particle representing the MBH was added at the cen- 
ter of each galaxy with a mass according to the updated 
Mbh~ A^buige relation of lHaring & Rixl (|2004|) . Our choice for 
the galaxy masses in conjunction with the assumption that the 
MBHs follow the MeH-A^buige relation result in MBH masses 
of 6 x 10 4 M and 6 x 10 5 M for the satellite and primary 
MBH respectively, with an initial MBH mass ratio q = q g!l \- 
With these choices, we target the typical masses and cosmic 
epoch of coalescing MBHs that should b e detectable by LISA 
dVolonteri et al.ll2003tlSesana et al.ll2005l) . 

We choose merger orbital parameters that are common fo r 
merging halos in cosmological simulations (Benson 2005). 
All simulations start with the baricenters of the two galaxies 
at a distance equal to the sum of their virial radii. Their or- 
bit is initially parabolic, with pericentric distance R p = 8 kpc, 
equal to 20% of the primary galaxy's virial radius, except for 
one case with a 10 times smaller pericenter. We also vary 
gas fractions and orbital inclinations. A summary of the sim- 
ulations here presented can be found in Table Q] 8 All sim- 

8 We note that only prograde mergers have been considered in this work. 



TABLE 1 

Summary of Merger Simulations 



f g primary 


f g satellite 


e b 


notes 


30% 


30% 


0° 


reference run 


30% 


50% 


0° 


1.7 X 10 5 SPH particles in the satellite a 


10% 


10% 


0° 




30% 


30% 


20" 




30% 


30% 


45" 




30% 


30% 


0° 


small pericenter (0.8 kpc) 



a This number has been chosen in order to have the same ratio of gas particle 
mass to MBH mass in both primary and satellite galaxies, as happens for all 
other simulations here discussed. 

b The orbital plane and satellite disk are inclined by an angle 8 with respect to 
the disk of the primary galaxy. 



ulations were performed w ith GASOLINE, a m ulti-stepping, 
parallel N-body /SPH code (IWadslev et al.ll2004l) . We include 
atomic cooling for a primordial mixture of hydrogen and he- 
lium, and the star formatio n algorithm is based on the lo- 
cal Schmidt-Kennicutt law dKatzlll992l) . Feedback from su- 
pernovae is t r eated using the blast- wave model described in 
IStinson et alJ d2006l) . Accretion on to the MBHs is also mod- 
elled with a sub-grid recipe (e.g.. ISpringel et al]|2005l) : the 
accretion rate Mbh is estimated from the density p g , sound 
speed c.5 and relative velocity V of the gas within a smooth- 
ing length from the MBH, via a Bondi-Hoyle-Lyttleton type 
formula, M BH = 4 7 rG 2 M 2 iH p i; (V 2 + c 2 r 3/2 dBondilll952l) . Of 
this mass-energy input, a fraction e r = 0.1 is assumed to be 
radiated away, while a fraction (1 — e r ) = 0.9 is added to the 
mass of the MBH from its neighboring gas particles. Lastly, 
a fraction eg, = 0.005 of the radiated luminosity couples to the 
surrounding gas as a heating source. This feedback efficiency 
is tuned so that a number of constraints re lated to our initial 
galaxy models can be satisfied (see Section [3~T1 >. In principle, 
the feedback efficiency depends both on the employed sub- 
grid model for the interstellar medium and the numerical res- 
olution. For all these reasons, our specific choi ce differs from 
those used in a number of earlier studies (e.g.JS pringel et ail 
120051) . Since the simulations are halted at satellite disruption, 
the orbital decay of the satellite MBH is determined by the 
well-resolved drag forces acting on its host galaxy; therefore, 
we do not include a sub- resolution Bondi dra g term, as em- 
ployed recently by, e.g.. lYounger et all d2008l) . Analogously, 
we do not employ any other numerical recipe that forces MBH 
orbital decay, because our aim is to study the growth of MBHs 
in relation to the efficiency of close MBH pair formation. In 
particular, the black holes were not repositioned at the local 
minimum of the grav itational potential at each time-step (e.g., 
Uohansson et 10120091) . 

3. RESULTS 

In Section lXTI our reference simulation is discussed, where 
the gas fraction in both dis ks o f th e pr imary and satellite 
galaxy is f s = 0.3. Sections 13.21 and [331 compare this refer- 
ence case with other merger simulations of different gas frac- 
tion and orbital parameters, respectively, as summarized in 
Table Q] 

3.1. Reference Simulation: Coplanar Merger with f g = 0.3 

In fact, retrograde encounters are characterized by a much slower orbital de- 
cay and, particularly in the low mass ratio regime, produce naked "wander- 
ing" MBHs as the satellite galaxy is disrupted at distances ~ 10 kpc. From 
such distances the dynamical friction timescale for producing a MBH pair at 
~ 100 pc are longer than a Hubble time. 
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FIG. 1. — Coplanar merger with a disk gas fraction of f g = 0.3 in both 
galaxies. Upper panel: Evolution of the mass of the MBH, Mbh, in the 
primary (dotted line) and satellite galaxy (solid line) as a function of time. 
Middle panel: Evolution of the mass ratio q of the two MBHs as a function 
of time. Lower panel: Orbital decay of the two MBHs as a function of time. 

In order to assess the effects of the merger on the mass 
growth of the MBHs, we first performed simulations of the 
primary and satellite galaxies evolved in isolation. For these 
tests, we chose the / g = 0.3 galaxy models of our reference 
simulation. In these simulations, the mass ratio of the two 
MBHs in the isolated galaxies does not deviate significantly 
or systematically from the initial q = 0.1 (the maximum fluc- 
tuation around this value is ~ 10%) over a period of more than 
2 Gyr. Therefore, we satisfy our working hypothesis that, in 
equilibrium conditions, the q of the black holes corresponds 
to the galaxy mass ratio, following from our initial choice 
of MaH-A^buige and galaxy morphology. In addition, these 
tests provide a measure of the MBH "quiescent" accretion: 
mass-doubling timescales of the MBHs evolved in isolation 
are ~ 2 Gyr, i.e. comparable to the typical duration of one 
of the mergers presented below. This result indicates that any 
larger MBH growth or large variation of q during the merger 
does not stem from secular evolution in the galaxy models or 
from numerical effects, but rather should be attributed to the 
galaxy encounter itself. 

Figure [Upresents the mass evolution of the two MBHs, the 
evolution of their mass ratio, and their relative separation as a 
function of time. By the end of the merger (f ~ 2.6 Gyr), ow- 
ing to dynamical friction, the two MBHs have formed a close 
pair in the nucleus of the remnant at a separation compara- 
ble to the adopted force resolution. This finding confirms our 
previous results (C09) and suggests that gas accretion onto 
the MBHs and associated feedback is not critical for pair for- 
mation in this case. The primary MBH grows quiescently 
throughout most of the merger, while the secondary increases 
its mass tenfold by the time the pair forms. The correspond- 
ing increase in the mass ratio of the two MBHs, q, occurs in 
two distinct episodes, which are elucidated in Figure[2] 

The upper panel of Figure [2] shows the evolution of the gas 
specific angular momentum in the direction of the disk rota- 




t [Gyr] 

FIG. 2. — Evolution of properties related to the coplanar merger with f g = 
0.3 in both galaxies. Upper panel: Gas specific angular momentum in the 
direction of the disk rotation axis in the primary (dotted line) and companion 
galaxy (solid line) (see text for details). Middle panel: Orbital eccentricity 
of the satellite MBH inside the disk of the primary. Lower panel: Eddington 
ratios, /kdd, of the mass accretion rates onto the primary (dotted line) and 
secondary MBH (solid line). 

tion axis, L z .g as /M gas , as a function of time for both galaxies. 
For this calculation, the gas particles that contribute to the 
baryonic mass within the central 5 softening lengths of each 
galaxy right before the third pericentric passage (f = 1 .7 Gyr) 
were traced back in time. This figure shows that most of 
the gas in the nuclear region of the satellite at this stage 
has lost its specific angular momentum on a relatively short 
timescale. The reason for the angular momentum loss is the 
strong tidal torques that occur near the second pericentric pas- 
sage (f = 1.2 Gyr) and are induced by the gravitational inter- 
action with the primary galaxy. As a result, a ~ 0.5 Gyr long 
accretion episode with a corresponding increase of q is ob- 
served. As shown in the lower panel of the same Figure, the 
accretion rate M onto the satellite MBH is 10% of the Edding- 
ton limit Mndd (Eddingt on! 19 161) during this phase. 

On the other hand, the gas that is funnelled near the center 
of the primary galaxy has been experiencing a nearly steady 
loss of specific angular momentum over a very long timescale 
(Figure |2j . This indicates that angular momentum loss in the 
case of the primary galaxy is not caused by tidal torques aris- 
ing from the interaction with the satellite, which would oc- 
cur at pericentric passages and become stronger as the merger 
progresses. Rather, it is induced by secular evolution (i.e. spi- 
ral arms) which redistributes angular momentum throughout 
the disk. In this context, the effect of initial transient spiral 
arms is evident during the first ~ 200 Myr. This agrees with 
the fact that mass growth of the primary MBH is essentially 
unchanged between the merger and the evolution in isolation. 

Around the third pericentric passage, ram pressure exerted 
by the interstellar medium of the primary galaxy strips all the 
gas from the satellite down to our force resolution. This is 
in agreement with an alytic estimates based on the study by 
iMarcolini et al.l d2003l) . as discussed also in C09. As a result, 
the satellite is now devoid of gas (as shown in left panel of 
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FIG. 3. — Gas density maps at t = 2.33 (left panel) and t = 2.48 Gyr (right panel), which correspond to just before and after the orbit of the secondary MBH 
circularizes inside the disk of the primary. The time span between these two snapshots corresponds roughly to two orbits of the satellite. The maps show the inner 
5 kpc of the primary. Density is projected onto the x—y plane and is color-coded on a logarithmic scale with brighter colors corresponding to regions of higher 
gas density. The black dashed line marks the trajectory of the satellite MBH and the square indicates the region around it. The satellite does not appear in the left 
panel, as ram pressure has stripped its entire gas content. On the contrary, the position of the satellite in the right panel is evidently traced by the overdensity and 
wake excited in the primary disk. 



Figure O, and accretion onto the smaller MBH is suddenly 
halted. A period of slowly decreasing q follows. In fact, dur- 
ing this phase, the more massive MBH continues to accrete 
gas from its host, experiencing an increase in its Eddington 
ratio / E dd = M BH /M E( jd as the satellite galaxy is now orbit- 
ing close enough to excite gas inflows in the primary disk 
(r>2Gyr). 

The mass of the satellite MBH sharply increases again (with 
an associated second increase in q) at kiloparsec-scale sepa- 
rations. At this stage, the secondary MBH orbits inside the 
gaseous disk of the primary. The mass increase coincides with 
a sudden drop in the orbital eccentricity of the satellite MBH, 
as shown in the middle panel of Figure|2] This drop in eccen- 
tricity is caused by dynamical friction acting on the satellite 
along its prograde coplanar orbit in the high-density region 
of the primary disk. Such orbit circularization is analogous 
to tha t found for MBHs in circumnuclear disks dDotti et aT] 
2009). Thus, the satellite MBH and its host stellar cusp are 
moving with a low relative velocity with respect to the disk 
of the primary. As a consequence, they are able to collect 
surrounding gas with low angular momentum in the reference 
frame of the satellite, creating an overdensity (right panel of 
Figure [3]) from which material is efficiently accreted by the 
satellite MBH up to a peak fm& ^0.3. On the other hand, ac- 
cretion onto the primary MBH still relies on angular momen- 
tum transport by instabilities in the disk. Such instabilities, 
triggered by the sinking satellite, are stronger than at earlier 
times, but still not able to sustain high /add- 

Overall, by the time the satellite is disrupted reaching the 
nuclear region of the merger remnant and the two MBHs form 
a close pair, the combination of strong tidal torques and orbital 
circularization acting on the companion galaxy has caused the 



MBH mass ratio to increase from 1 : 10 to 1 : 3, bringing it into 
a regime of "major" mass ratio. 

3.2. The Role of Gas Fraction 

Figure [4] compares the evolution of the mass ratio of the 
MBHs q as a function of time in mergers with different disk 
gas fractions, f g (see Tabl^T)). All of these mergers are on 
prograde, coplanar orbits. Again, we follow the evolution of 
the interacting systems up to the time when the orbital decay 
of the satellite galaxy is completed. This figure shows that 
the first phase of increase in q, caused by dynamical destabi- 
lization of the satellite, happens in all three cases: due to the 
fact that inclination and orbit are fixed in these mergers, the 
torques acting on the satellites have the same strength. Inter- 
estingly, the highest value of q reached in the first stage traces 
roughly the amount of gas oc / g available for MBH fuelling. 
In all the mergers considered here, ram pressure is effective 
in removing gas from the satellite galaxy, once the two galaxy 
disks come into contact. When this happens, accretion onto 
the secondary MBH is halted. 

Figure|4]shows that the initial / g does cause significant dif- 
ferences in the evolution of q also during the second stage 
of the mergers. Indeed, the second phase of strong accre- 
tion onto the satellite MBH is absent in the merger with the 
smallest gas fraction (f g = 0.1). As discussed earlier, mass 
growth during this phase becomes efficient only when the sec- 
ondary MBH moves inside the gas disk of the primary galaxy 
and its orbit circularizes. Instead, in the case of f g = 0.1 the 
satellite does not sink below ~ 400 pc before being tidally 
disrupted, because it has not experienced strong gas inflows 
steepening its potential well during the first orbits (see C09 
for a detailed discussion). At such large distances, the MBH 
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FIG. 4. — Evolution of the mass ratio q of the two MBHs as a function of 
time in coplanar mergers with different gas fractions / g . Dotted (blue) and 
solid (red) lines show results for the mergers with gas fractions fg =0.1 and 
fg = 0.3, respectively, in both galaxies. The dashed (green) line corresponds 
to a merger where the initial gas fractions are f g = 0.3 in the primary and 
fg = 0.5 in the satellite galaxy. 



orbit is still mildly eccentric, and the background density is 
not high enough to allow for the formation of an overdensity 
that could trigger the second accretion episode. As a conse- 
quence, the MBH mass ratio at the time of satellite disruption 
is 1 :6, close to its initial value, and the naked satellite MBH 
will still take a few billion years to sink to the center of the 
remnant via dynamical friction (C09). 

We now focus on the merger where the initial gas fraction 
in the primary and companion galaxies is equal to / g = 0.3 
and / g = 0.5, respectively. Figure |4] shows that this case is 
characterized by a much larger final increase in q compared 
to our reference case, where the initial gas fraction was equal 
to / g = 0.3 in both galaxies. Bearing in mind that the only dif- 
ference between the two initial conditions is the satellite / g , 
and that the satellite gas has already been entirely stripped by 
ram pressure at this late stage, this interesting result can be 
explained as an effect of the different structure of the surviv- 
ing satellite core and the different mass growth of its MBH. In 
fact, a larger initial gas fraction allows the satellite to build a 
denser stellar core via star formation in response to tidal per- 
turbations during the first two orbits. Consequently, the nu- 
clear region of the satellite is denser and more massive. It is 
therefore more efficient at collecting gas from the disk of the 
primary and it is subject to an enhanced sinking. As a result, 
the orbit of the secondary MBH undergoes circularization in a 
denser region of the primary disk, and is able to accrete more 
gas. Moreover, by the time the secondary MBH enters the 
disk of the primary, it is ~ 60% more massive compared to the 
case of f s = 0.3. This difference in mass naturally enhances 
its accretion rate, which scales as oc M| H . At the end of this 
merger, the satellite MBH has sunk down to the center of the 
remnant, where a close MBH pair is formed at scales ~ 30 pc, 
comparable to our force resolution. The mass ratio of the pair 
at this point is the highest among all our simulations: q = 0.5. 
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FIG. 5 . — The evolution of orbital parameters as a function of time is shown 
for the inclined runs: the continuous (red) line shows data for the merger with 
20° inclination, the dashed (blue) line for the one at 45° . Upper panel: MBH 
separation. As expected, sinking is more effective for more coplanar orbits. 
Middle panel: angle 9 between the orbital plane of the two MBHs and the 
plane of the primary galaxy's disk. Orbit dragging onto the disk plane can 
be seen in both cases, once the satellite's orbit reaches the disk-dominated 
region of the primary galaxy. Lower panel: Mass ratio of the MBHs. 

3.3. The Role of Orbital Parameters 



As shown in Sections 13.11 and 13.21 the evolution of the 
MBH mass ratio depends on processes (especially ram pres- 
sure stripping of gas and orbit circularization) that may be 
sensitive to the initial orbit of the merger itself. In this Sec- 
tion, we discuss results from three simulations addressing the 
dependence on orbital parameters. In all of them, the gas frac- 
tion in the two galaxies was equal to that of the reference case, 
/„ = 0.3. 

Upper panel of Figure [5] shows the evolution of the MBH 
orbital separation in two merger simulations where the or- 
bital plane and satellite disk were inclined with respect to 
the plane of the primary disk, by 20° in one case, and 
45° in the other. Torques acting on the satellite during the 
early phases of the merger are weaker for higher inclinations 
(Barne s~& Hernquistlll996l) . and for this reason the increase 
in mass ratio q during the first three orbits is milder than in 
the coplanar case with the same gas fraction (Figure|5] lower 
panel). Moreover, a higher inclination corresponds to a slower 
orbital decay; in fact, the satellite spends only a small fraction 
of each orbit in the high density, co-rotating baryonic region 
of the primary galaxy, where dynamical friction is most ef- 
ficient. The orbi t is eventually dragged do wn to the plane of 
the primary disk (lOuinn & Goodman! 19861) (see middle panel 
Figure |5]). However, since such drag takes a number of or- 
bits to cause significant alignment, the circularization effect 
- which can only happen when orbiting for most of the time 
inside the primary disk - is delayed compared to the copla- 
nar case. For this reason, the satellite undergoes a larger 
number of tidal shocks and is eventually disrupted, preclud- 
ing a second episode of substantial accretion onto the satel- 
lite MBH. For low initial inclination (9 = 20°), the drag is 
effective enough to bring the MBH down to a separation of 
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FIG. 6. — MBH evolution in the fg = 0.3 merger with 0° inclination and 
small pericenter. Upper panel: mass accretion of the MBHs. The growth of 
the satellite MBH stops at second passage (f = 0.8 Gyr), earlier than in the 
other cases, as ram pressure stripping of the satellite gas is more efficient at 
smaller separations. The primary MBH, on the other hand, starts a period of 
enhanced accretion after the second orbit, as the perturbations excited onto 
the primary galaxy are correspondingly stronger. Lower panel: evolution of 
the MBH mass ratio q. 

~ 120 pc before the satellite is disrupted. If the inclination is 
higher (8 = 45°), the satellite is disrupted leaving the MBH at 
~ 700 pc distance, where the dynamical friction timescale for 
sinking of a naked MBH is of a few billion years. 

Finally, we turn to the merger starting on a parabolic orbit 
with very small pericenter (0.8 kpc). The MBH masses and 
q as a function of time in this case are shown in Figure [6] 
During the first orbits, the satellite in this merger is subject to 
stronger tidal torques than in the reference case, but MBH ac- 
cretion is weaker. This is because the ram pressure acting on 
the satellite's gas is P oc p ex tV 2 , where p ext is the gas density 
in the primary disk and V is the relative velocity between the 
satellite and its surroundings. Because of the small distance 
(corresponding to a high p ext ) and high V at pericenters, ram 
pressure stripping is effective since the very first passage: the 
satellite MBH is starved early and is unable to accrete a sub- 
stantial amount of mass. For the same reason, MBH growth in 
the satellite is halted completely already at second pericenter 
(t ~ 1.8 Gyr). On the contrary, the primary MBH experiences 
enhanced accretion: in fact, the close passages of the satel- 
lite excite strong instabilities in the primary disk, funnelling 
gas towards its nucleus. For this reason, at intermediate times 
during the merger, the mass ratio drops down to 1 : 20. Anal- 
ogously to the reference run, a second phase of accretion onto 
the satellite MBH occurs, in this case bringing q back to its 
starting value. However, MBH sinking does not proceed any 
further: as a consequence of the strong tidal shocks experi- 
enced by the satellite at pericentric passages, the satellite is 
disrupted leaving its MBH at a distance of 1 kpc. At these 
separations, the estimated timescale for orbital sinking of the 
naked MBH is longer than a Hubble time. 

4. DISCUSSION AND CONCLUSIONS 
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FIG. 7. — The different symbols show the separation of the two MBHs 
and the corresponding MBH mass ratio q at the time of satellite disruption 
for all the 1:10 mergers discussed here, labelled according to their initial gas 
fractions f g , orbital inclination 8 and initial pericenter R p . 

The outcomes of the simulations presented in this Paper are 
summarized in Figure|7l where the mass ratio q of the MBHs 
is plotted against their separation, at the time of satellite dis- 
ruption. As discussed in the previous Sections, the physical 
processes that facilitate the MBH orbital decay correspond to 
a stronger accumulation of gaseous mass in the central re- 
gion of the satellite, compared to the isolated, equilibrium 
galaxy models. In turn, such gas density enhancement cor- 
responds, in our modelling of black hole accretion, to an in- 
crease in the MBH mass ratio. Applying the s tandard dynam- 
ical friction formulation ( Chandrasekhari n"943h to a "naked" 
MBH in our merger remnants, we estimate the timescale for 
bringing the MBH down to the nuclear region (< 30 pc) to 
be up to a few billion years for separations up to ~ 1 kpc, 
and comparable to or longer than a Hubble time for larger 
separ ations. Thanks to a combina tion of stellar-dynamical 
(e.g..lMilosaylievic & Me rritt 2001) a nd gas-dynamica l (e.g., 
iGould & Rixll2000HEscala et alj|2005t iDotti et al.ll2007l) pro- 
cesses, a satellite MBH that sinks down to our current force 
resolution limit may decay further and form a Keplerian bi- 
nary together with the primary MBH. Therefore, our results 
suggest that mergers that may most promptly produce MBH 
binaries are also those that tend to enhance the MBHs q, with 
respect to what inferred from the galaxy mass ratio q ga i. 

Such increase in MBH q happens in two distinct phases, 
whose occurrence and relative importance depend on the de- 
tails of the merger process itself. Specifically, in the initial 
stages of the encounter, the stronger tidal perturbations ex- 
perienced by the satellite galaxy, compared to those of the 
primary, cause an enhanced mass growth of its MBH. In ad- 
dition, in the last stages of the encounter, the orbit of the sec- 
ondary MBH may circularize inside the disk of the primary. 
As a result, its ability to accrete gas and grow in mass relative 
to that of the primary MBH can be further amplified. Such cir- 
cularization and associated increase in the accretion rate has 
been previously reported in small-scale simulations of MBH 
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pairs embedded in a common nuclear disk dDotti et al.ll2009l) . 
We note that the amount of gas left around the satellite MBH 
by ram pressure stripping might be underestimated. In fact, 
our simulations do not resolve the accretion disk around the 
MBHs, nor the cold molecular phase in the nuclear region of 
the satellite; both of th ese would be less susceptible to ram 
pressure stripping (e.g., Quilis et al. 2000) and could remain 
bound to the satellite while the rest of its ISM is lost. In this 
case, the growth of the secondary MBH may not stop com- 
pletely and the mass ratios may become even larger. 

A cautionary remark concerns the fact that MBH accre- 
tion has been modelled assuming the Bondi-Hoyle-Lyttleton 
sub-grid accretion recipe that is widely used in the literature. 
However, accretion will most likely occur by means of an- 
gular momentum tra nsport in a nuclear d isk with an effec- 
tive a-viscosity (e.g.. lLin & Pr ingle 1987). Other models for 
gas funnelling and accretion onto central MBHs have been re- 
cently proposed to account for the shortcomings of the "stan- 
dard" approach. iPower et"aT1 (1201 Oh propose an accretion 
scheme for angular momentum-dominated gaseous flows in 
a stationary gravitational background, introducing a delay in 
the accretion as a free parameter related to viscous timescales; 
from their simulations, they find that the Bondi prescription 
overe stimates the inflow rate acr oss the MBH influence ra- 
dius. iHopkins & Ouataertl d2010l) model the transport of an- 
gular momentum via disk instabilities arising in galaxy merg- 
ers, on scales both larger and smaller than those resolved in 
this work. They describe the inflow rates down to the MBH 
sphere of influence with M ~ aM g £l, where M g is the avail- 
able gaseous mass, £1 the angular speed, and a is related to 
the amplitude of the instabilities. They find that the Bondi ac- 
cretion recipe systematically underestimates (although with 
large scatter) the actual inflow rate at sub-pc scales, while 
better predictions are given by an effective a-disk model. It 
is therefore unclear whether our accretion recipe introduces 
a bias in the MBH mass growth, and in which direction. In 
an attempt to investigate the effect of other sub-grid recipes 
on our results, we computed a posteriori accretion rates in 
our refere nce simulation with th e recent a-disk sub-resolution 
model of iDeBuhr et ail d2009l) . We find that such model 
yields higher M compare d to the Bondi prescription, in agree- 
ment with the findings by Hopkins & Ouataertl (1201 Ol) ; more- 
over, the accretion rates onto the two MBHs are enhanced by 
roughly the same factor. While the absolute values of M can 
depend on the employed numerical prescription, the physi- 
cal picture concerning the relative growth of the two MBHs 
emerging from our simulations may not be substantially af- 
fected by the specific choice of sub-grid modelling. Indeed, 
our findings reflect clear and well-resolved large scale effects, 
namely how gravitational torques and orbit circularization can 
enhance MBH sinking while making a larger gaseous mass 
relative to Mbh available for accretion onto the secondary 
black hole. 

The results presented in this paper are especially 
relevant in the context o f MBH gravitational recoils 
(lLousto & Zlochowerl l2009t iTanaka & Haimanl 120091) . In- 
deed, if a large fraction of unequal-mass galaxy mergers re- 
sults in mergers between MBHs with nearly equal masses, 
then the recoil velocity distribution o f the MBH populatio n 
will be different than expected (e.g.. IVolonteri et aD l2010h . 
However, the actual recoil velocity distribution will also de- 
pend on the magnitude and relative orientation of the spins 
of the MBHs at the final stage of the merger, which is likely 



driven by gas dy namics at scales well below those res olved in 
our simulations dPerego et al.ll2009b iDotti et alJl2010h . 

Our findings, togeth er with those of C09 and of 
iKazantzidis et al.l d2005l) . confirm the fundamental role of 
tidal stripping and gas-dynamical effect in deciding the for- 
mation and properties of close MBH pairs in unequal-mass 
galaxy mergers. Moreover, our results suggest that the ef- 
ficiency of MBH pair formation may correlate with the fi- 
nal mass ratio of the pair itself, so that MBH pairs with 
larger mass ratios tend to be produced more effectively and 
promptly. Gravitational wave detectors, such as LISA, will 
enable the use of gravitational wave signals from MBH coa- 
lescences as a new, independent probe of cosmic structure for- 
mation. Indeed, gravitational wave-forms can allow the deter- 
minatio n of mass, spi n, and orbital parameters of the merging 
MBHs ( Vecchio 2004). In principle, this information could be 
used to infer the masses of the merging host galaxies. How- 
ever, our results show that this connection cannot be made by 
simply applying the observed scaling relations between the 
masses of MBHs and the properties of their host galaxies, 
even if these scalings are applicable to galaxies throughout 
cosmic history. The findings presented here suggest that the 
mapping between galaxy and MBH mergers depends on var- 
ious factors, such as the gas content of the merging galaxies 
and the encounter geometry, and as such might need to be ap- 
proached in a probabilistic way. Such investigations would 
require a combination of a series of merger experiments that 
explore a larger parameter space with semi-analytical models 
of the co-evolution between galaxies and MBHs. 
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